Effective Mass of Atom and Excitation Spectrum in Liquid 

Helium-4 at T = K 

A. A. Rovenchak 

Department for Theoretical Physics, Ivan Franko National University of Lviv 
12 Draghomanov St., Lviv, UA-79005, Ukraine 
tel: +380 322 979443; e-mail: andrij@ktf.franko.lviv.ua 

February 1, 2008 

Abstract 

A self-consistent approach is applied for the calculations within the two-time temperature 
Green functions formalism in the random phase approximation. The effective mass of 4 He 
atom is computed as m* = 1.58 m. The excitation spectrum is found to be in a satisfactory 
agreement with the experiment. The sound velocity is calculated as 230 m/s. The tempera- 
ture of the Bose-condensation with the effective mass taken into consideration is estimated 
as 1.99 K. 
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1 Introduction 

The idea of the effective mass of helium atom was suggested by Feynman in 1953 [|TJ. He stated 
that one should insert the effective mass (being slightly larger than the mass of a 'pure' atom m) 
in the expressions for the density matrix. 

Isihara and Samulski used the effective mass m* = 1.71m in order to obtain a good 
agreement of the excitation spectrum sound branch with the experimental data on the sound 
velocity. 

While in both these cases the effective mass was introduced phenomenologically, it appeared to 
be possible to receive the value of this quantity basing on the experimental data for the structure 
factor of liquid helium-4 ||, the result of 1.70 m was calculated there by Vakarchuk. 

In the paper presented we will give the way of the receiving 4 He atom effective mass by means of 
a self-consistent equation. The expressions are written within the collective variables formalism as 
described by Bogoliubov and Zubarev H. Two-time temperature Green functions || are utilized 
for the calculation of the thermodynamic averages. 

The main idea of the paper is to show the possibility of an essentially simple approach to the 
problem of the many-boson system with strong interaction, such as liquid helium-4. The method 
applied does not require much computational efforts. This advantage allows for the development 
of the further approximations. 

In addition, if one accepts the assumption that the phenomena in liquid 4 He are at least 
partly due to the Bose-condensation being 'spoiled' by the interatomic interaction, it turns to be 
possible to estimate the lambda transition temperature as the critical temperature of the ideal 



Bose-gas. We show that such an approach leads to a very good agreement with the experiment: 
if the effective mass is about 50% larger than the pure one the Bose-condensation temperature 
decreases to the value of ~ 2 K. 

The Green functions technique provides also a possibility to receive the excitation spectrum 
of the system. As a result of the mass renormalization, the excitation spectrum is found to be in 
the better agreement with the experiment in comparison with the Bogoliubov's or Feynman's one 
while the expressions are the same (for the latter two spectra the problem of the so called 'roton' 
minimum overestimation is well-known if one considers the pure mass). 

A self-consistent approach was recently applied for the calculation of the 4 He excitation spec- 
trum by Pashitskii et al 0. The results of this paper are in an excellent agreement with the 
experiment. The authors used the 'semitransparent spheres' potential for the calculation with 
some adjusting parameters. In our work, we utilize the potential as an input information for 
the computations. This potential was received on the basis of the quantum-mechanical equations 
with the static structure factor as the only experimental data. Since this quantity is quite easily 
measured directly in the scattering experiments, we consider it as a good approach. Unfortu- 
nately, a direct calculation of the potential for the many-body problem cannot be made for the 
time being. 

The paper is organized as follows. The calculating procedure is given in Section 2. The 
Hamiltonian is written and the equations of motion for the Green functions are solved in random 
phase approximation (RPA) providing a self-consistent equation for the effective mass extraction. 
The numerical results are adduced in Section 3 together with the discussion. 



2 Calculation procedure 

The Hamiltonian of the Bose-system in the collective variables representation reads ||J]: 
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where the operator <9k = <9/<9p_k- Here, Sk is the energy spectrum of a free particle, Sk = 
h 2 k 2 /2m, N is the total number of particles in the system, and V is the system volume. In the 
thermodynamic limit, N/V = g = const. Vk is the Fourier transform of the interatomic potential. 
The item with one summation over the wave vector k in Eq. ([I]) corresponds to the random phase 
approximation, and the second one is the correction. Let us assume that our system is described 
by exactly RPA Hamiltonian H^*\ i. e., 
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where el = h 2 k 2 /2m* and m* is the effective mass of 4 He atom. It is the only quantity suitable for 
the 'effective' role since we wish to preserve the interatomic potential as the initial information. 

One can define m* demanding that the effective Hamiltonian (0) leads to the same ground-state 
energy as the initial Hamiltonian ([!]), (H^) = (H): 
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where the superscript (*) near the angle brackets is introduced for the convenience. 

One can find the operators product average utilizing two-time temperature Green functions 
defined as follows 0: 

((A(t)\B(t'))) = i9(t-t')([A(t),B(t')]) (4) 
with operators given in the Heisenberg representation, 6 is the Heaviside step function. 
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where Gba stands for ((B\A)) and the operator A is introduced for the convenience. We put 
the time arguments in the operators A(t),B(t') to coincide: t — t' = 0. This will provide static 
properties of the system under consideration. In the above expression, j3 is the inverse temperature, 
P=l/T. 

Now, we will proceed to the equations of motion for the Green functions G pp (k) = ((pk|p-k))i 
G p a(k) = ((pk|<9-k)), etc. It is easy to receive the following set of equations in RPA: 

(hu + e k )G pp (k) = 2e k G dp (k), 

(hu - e k )G dp (k) = QVkG pp (k) + — , 

(hu + e k )G pd (k) = 2e k G dd (k) - —, 
(Ok) - e k )G dd (k) = Qv k G pd (k). 
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Next, let us consider the triple product average (ABC). One can obtain it utilizing either the 
Green function Gc-ab = ((C\AB)) or Gbc-,a = {{BC\A)). We suggest the first possibility to fulfill: 
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In other words, we neglect the functions of the type Gbc,a = ({BC\A}) for the sake of simplicity 
(when applying this to Eq. (D) it means that only the RPA term of the Hamiltonian (|l|) is taken 
into consideration when constructing the equations of motion). 

Having performed the similar procedure with the function Gg ;p g(ki, k 2 , k 3 ) = ((d\ il \p\ i2 d\ ! _ 3 )) we 
obtain in RPA the following set of equations: 
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where the quadruple Green functions were decoupled in such a way providing for the inhomoge- 
neous set of equations: 



((AB\CD)) = (BD)((A\C)) + (CA)((B\D)) + (AD)((B\C)) + (CB)((A\D)). 
The inhomogeneous terms in Eq. (§) read: 
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The notations for the averages of pair products are listed below: 
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The quantity a k is defined as follows: 
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Now, if we turn back to correlation (Q), the meaning of the asterisk as a superscript becomes 
clear: one should substitute m with m* in the l.h.s. of this equation. 

In the ground state (T = K), hyperbolic cotangents in Eq. ( |TTD equal to unity Therefore, a 
self-consistent equation for the extraction of m* becomes as follows: 
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where p = k + q. We also consider the specific energy instead of the total one by introducing the 
factor of 1/JV. 



3 Numerical results and discussion 

We use the previously obtained results for the interatomic potential Fourier transform v k . 
The value of the equilibrium density is g = 0.02185 A" 3 . The mass of helium-4 atom equals 
m = 4.0026 a.m.u. We transit from the summation over the wave vector to the integration in the 
usual way: ^ k — » V J dk/(27r) 3 . The value of the upper cut-off for the integration over the wave 
vector is 16.0 A -1 . 

The solution of Eq. (|T3| ) at the abovelisted conditions is 

m* = 1.58 m. (14) 

One can also receive the excitations spectrum using Green functions. The solutions of set (^) 
are proportional to l/(ti 2 ui 2 — £ k o> k ) providing the spectrum E k = ±e k a k — a very well known 
result [Q. If one inserts the effective mass into the definitions of e k and a k the obtained curve fits 
the experimental one quite satisfactory, see Fig. [I]. 



The phonon branch is reflected quite good providing the sound velocity of approximately 
230 m/s vs the experimental 238 m/s at T = 0.8 K @ or 240 m/s at T = 0.1 K 0. The so called 
'roton' minimum also has the value close to the experimental one. 

In addition, the obtained value of the effective mass shifts the temperature of the Bose- 
condensation from T c = 3.14 K for the pure mass to T c = 1.99 K vs the experimental temperature 
of the lambda transition T A = 2.17 K. We consider the results discussed above as quite good as 
for such a rough approximation as random phases. 




Figure 1: Excitations spectrum of liquid helium-4. Filled circles — experimental data ||; Solid 
line — calculated energy. 
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